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A relativistic constituent quark model is found to reproduce the recent 
data regarding the ratio of proton form factors, F2{Q'^)/Fi{Q'^). We show 
I that imposing Poincare invariance leads to substantial violation of the helicity 

conservation rule, as well as an analytic result that the ratio F2{Q'^) /Fi {Q"^) ~ 
1/Q for intermediate values of . 
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> \ I. INTRODUCTION 

(N 

O I The recent discovery [|l|, [0] that the ratio of proton form factors Ge/Gm falls linearly 

^ ' with and that the ratio QF2/F1 reaches a constant value for 2 < < 6 GeV^ has 
CN ■ focused attention on understanding nucleon structure. The constant nature of the ratio 

^ ' QF2/ Fi contrasts with the prediction from perturbative QCD |^,^ that Q^F2/Fi should be 
constant. While this latter ratio could be achieved when experiments are pushed to higher 
values of Q^, it is worthwhile to obtain a deeper understanding of the present results. In 
particular, the perturbative result is based on the notion that helicity is conserved in 
Ch ■ high momentum transfer exclusive processes, so that it becomes interesting to understand 
why this conservation does not seem to be applicable. 
I The qualitative nature of the experimental results were anticipated or reproduced by 

I several model calculations , with that of Ref . ^ based on the work of Schlumpf |IO|-p!^ 
being the earliest. The implementation of relativity is an important feature of each of these 
calculations, so it is natural seek an understanding of the form factors in terms of relativity. 
Our purpose here is to examine the model of Ref. with the aim of highlighting the essential 
features which cause the ratio QF2/F1 to be constant. 

We proceed by presenting definitions and kinematics relevant for a light front analysis 
in Sect. II. The relevant features of our relativistic constituent quark model are displayed in 
Sect. III. The essential reason for the constant ratio QF2/F1 is discussed in Sect. IV, and 
elaborated upon numerically in Sect. V. The paper is concluded with a brief summary. 



II. DEFINITIONS AND KINEMATICS 

The electromagnetic current matrix element can be written in terms of two form factors 
taking into account current and parity conservation: 
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with momentum transfer = {p' —p)'^, = —q^ and is taken as the electromagnetic 
current of free quarks. For = Q the form factors Fi and HF2 are respectively equal to 
the charge and the anomalous magnetic moment k in units e and e/Mjv, and the magnetic 
moment is /i = -Fi(O) + kF2(0) = ! + «;. The Sachs form factors are defined as 



KF2 , and Gm = Fi + KF2 



(2) 



The evaluation of the form factors is simplified by using the so-called Drell Yan reference 
frame in which g"*" = so that = l\- This means that the plus components of the nucleon 
momenta (and also those of the struck constituent quark) are not changed by the absorption 
of the incoming photon. 

If light-front spinors for the nucleons are used, the form factors can be expressed simply 



in terms of the plus-component of the current [13]: 
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The form factors are calculated, using the "good" component of the current, J"*", to suppress 
the effects of quark-pair terms. 

It is worthwhile to compare the formalism embodied in Eq. @ with the non-relativistic 
quark model formalism in which Ge and Gm are the Fourier transforms of the ground state 
matrix elements of the quark charge (X]j=i,3ej5(r — ri) and magnetization I]i=i,3 "i^^i^ 
density operators. At high momentum transfer one needs to account for the influence of 
the motion of the proton on its wave function. Since the charged quarks of the initial 
proton state (or final state, or both initial and final states) are moving, charge and magnetic 
effects are correlated in a manner consistent with relativity. It is necessary to maintain 
this relativistic connection, lost in the non-relativistic quark model. The use of light front 
dynamics, concomitant with Eq. (|^), is a particularly convenient way to handle the motion 
of the initial and final state. This is because the proton wave function, a function of internal 
relative coordinates, is the same in any reference frame. 



III. RELATIVISTIC CONSTITUENT QUARK MODEL OF THE NUCLEON 



We study the form factors using relativistic constituent quark models in general, and 
starting with the model of Schlumpf p0|-p!2| in particular. Such models have a long his- 
tory [14-16|, and many authors 0-32| have contributed to the necessary developments. 
Schlumpf 's model is used because his power-law wave functions lead to a reasonably good 
description of the proton electromagnetic form factors, Ge and Gm-, at all of the values of 

where data were available as of 1992 [|TU]-[T^. This model uses the Bakamjian-Thomas 
BT construction, which implies the choice of a very specific model wave function 



23,31 



The use of a definite model allows us to gain insight, but we also shall discuss the limitations 
of this approach. 
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We remind the reader about a few basic features of light front treatments in the BT 
approach. The hght-front formahsm is specified by the invariant hypersurface = x^+x^ = 
constant. The following notation is used: A four-vector is given by A^^ = A^ , Aj_), 
where A^ = A^ ± A^ and A_|_ = {A^jA"^). Light-front momenta vectors are denoted by 
P = (P^?P±)) with p~ = + Tn?)/p^ for on-shell quarks. The three momenta pi of the 
quarks can be transformed to the total and relative momenta to facilitate the separation of 
the center of mass motion as 

o , , c Pt+Pt 

P = pi + P2 + P3, ^= —^—^ — + , V = — HT— ' 

Pi + P2 

k± = (1 - Opi± - CP2± , = (1 -r7)(pi_L + P2±) - ?7P3± • (4) 

It is also useful to consider the mass operator of a non-interacting system of total momentum 

where m is the light quark mass, taken as the same for up and down quarks. 

One may express the proton wave function in the center of mass frame in which the 
individual momenta are given by 

Pi± = k^ + eKx, P2± = -k^ + (1 - 0K± , P3± = -K^. (6) 

The use of light front variables enables one to separate the center of mass motion from the 
internal motion. The internal wave function \1/ is therefore a function of the relative momenta 
Pi, (^,7]. The internal wave function of the proton depends on these relative momenta. One 
can obtain the wave function in certain special frames via a kinematic boost. In particular, 
if the proton acquires a transverse momentum (is boosted) by the absorption of a photon 
of momentum q = (0, q^) by the third quark, the effects of the boost are obtained merely 
replacing the momenta k^, by 

kl = kx, = - m± (7) 

The mass operator, Mq of the boosted system is obtained by replacing k_L,K_L by the 
variables of Eq. so that 

^0 -M,+ ^^^^ . (8) 

We now turn to the construction of the non-perturbative wave function, This is based 
on the attempt to construct a state, described in terms of the given light front variables, that 
also is an eigenstate of angular momentum EO]. We take the proton wave function to be a 



product of an anti-symmetric color wave function with a symmetric flavor-spin-momentum 
wave function \I'. To understand the construction |^0| of the relativistic wave function, it is 



worthwhile to start by considering the non-relativistic model. Then 



^^« = i=(0,Xp + 0AXA)$ (9) 
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where 0p represents a mixed-antisymmetric and 0a a mixed-symmetric flavor wave function 
and, Xp,\ represents mixed symmetric or anti-symmetric spin wave functions (in terms of 
Pauli spinors). In the non-relativistic model the wave function $ depends on spatial variables 
only, and the computed form factors Ge and Gm will have the same dependence on Q'^. 
The relativistic generalization of Eq. (^) is 

■^ipi) = u{pi)u{p2)u{ps)ip{pi,p2,P3)y (10) 

where pi represents space, spin and isospin indices: pi = PiSi, Ti and repeated indices are 
summed over. The spinors u are canonical Dirac spinors: 

s) = i\ , (11) 

,jE{p) +m\ U / 

with the isospin label suppressed. 

The completely symmetric nature of the space-spin-flavor wave function is preserved by 
using 

Hpi,P2,P3) = $ [uip,)Tu^ip2) uMum{0)+u{p^)T^^''u^{p2) uMnM{0)g^J\ (12) 
r ^ --Li±^^5«r2, T'^'" ^ l=YCtT^T2, f^'-^^Y^^r^, (13) 

where the charge-conjugation matrix C = i'-f'^'-f^ = —i'~f5a2- Note that if one takes the non- 
relativistic limit of ([T2|) by taking pj — > (m, 0_l) one gets (|^) for the spin-isospin dependence. 

The momentum wave function can be chosen as a function of Mq to fulflU the require- 
ments of spherical and permutation symmetry. We take the S-state orbital function <I>(M) 
to be of a power law form: 

which depends on two free parameters, the constituent quark mass and the conflnement 
scale parameter /?. Schlumpf's parameters are /3=0.607 GeV, 7 = 3.5, and the constituent 
quark mass, m=0.267 GeV. 



The wave function of Eq. (jl^) now specifled. It is an eigenstate of the angular momentum 
operator (if interactions present in that operator are neglected) and also corresponds to a 
state of vanishing orbital angular momentum. This state is Poincare invariant according to 



Ref. pOl, but this invariance may be incomplete ||32 



IV. THE ESSENTIAL EFFECT 



The form factors of Eq. (|^) are obtained by computing the matrix elements of the current 
operator J"*". Here we take the quarks to be elementary particles, so that the operator J"*" 
is essentially the operator 7"*" times the charge of the quarks. The calculations may be 
simplifled by making a unitary transformation which replaces the Dirac spinors of Eq. (p!0D 
by light front spinors 
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(15) 



because 



«l(p+, p', X'h+ULip-^, p, A) = 26xyp+. 



(16) 



One then uses the completeness relation, 1 = J2x ul{p, ^)ul{p, A)/2m, in Eq. (|T0[), to obtain 
the light front representation for the wave function: 



^(Pi) = Ul{pi, Xi)ul{p2, X2)Ul{P3, X3)'lpL{Pi, Aj), 
i^L{Pi,Xi) = (Al|7^L(pl)|Sl)(A2|7^i/(p2)|s2)(A3|7^]^,/(p3)|s3)V^(pl,P2,P3), 



where TZm is a. Melosh rotation |^3[ acting between Pauli spinors. For example, 



(A3|7^]^,/(P3)|S3) = Ul(p3, A3)u(p3, S3) = (A3 



m + (1 - ?7)Mo + icr ■ (n X P3) 

m + {l-r])Moy+pl^ 



(17) 
(18) 



ks). (19) 



The net result of this is that the relativistic spin effect is to replace the Pauli spinors 
of Eq. (^ with Melosh rotation operators acting on the very same Pauli spinors [^]. The 
spin-wave function of the ith quark is given by 



T Pi) = ^M(Pi 



and I T Pi) = ^m(P 







(20) 



This means that, for example, the spin wave function xp is replaced by a momentum- 
dependent spin wave function Ix^f^Pi) 



XTP^) 



71' 



fPliP2 - iPltP2)|TP3) = 



lxr(Pl,P2))|TP3) 



(21) 



for a spin +1/2 proton. The term Xo^ represents the relativistic generalization spin-0 wave 
function of the quark pair labeled by (1,2). Note that the momenta Pi are to be expressed 
in terms of the relative variables k^, K^, k'j_, K'^ of Eqs.(§) and (|^). We shall see that the 
important relativistic effect is contained in the difference between the spinors of Eq. ( PO]) 
and Pauli spinors. 

The next step is to simplify the calculation by using the symmetry of the wave function 
under interchange of particle labels to replace the quark current operator by three times 
that operator acting only on the third quark. The average charge of the third quark of 
the mixed-symmetric flavor wave function vanishes. This means that the second term of 
Eq. ([1^ ) does not contribute to proton electromagnetic form factors. The first term involves 
a mixed-antisymmetric wave function, so the third quark carries the spin of the proton, S3 
of Eq. (|19]). However, the helicity of the light front spinor A3 can either be the same as or 
different from S3. The weighting of the terms is determined by the two terms of the Melosh 
transformation of Eq. (piOl) . Considering the arguments of the final state wave function allows 
one to understand that the two terms are comparable. The effect of the boost is incorporated 
simply by using Eq. (|^) in the argument of the final state wave function. Thus the large 
momentum Q is involved and both terms of the Melosh rotation (|T9|) are comparable. 
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There is clearly a substantial amplitude for a spin-up valence quark to carry a negative 
light front helicity. No suppression of light-front helicity-fiip can be expected from the 
use of such a wave function. This means that helicity conservation does not occur. This 
conclusion has been obtained by Ralston and collaborators based on the presence of 



non-zero orbital angular momentum. Here we have no orbital angular momentum, and the 
mixture of light-front helicity we obtain occurs as a result of imposing Poincare invariance 



on the constituent quark model. Nonetheless, we certainly support their statements 
that helicity non-conservation is an important effect. We also note that Braun et al. ||35 
argue that soft non-factorizable terms in the wave function, with helicity structure similar 
to ours, are important at intermediate values of . In our model, there is no basis for 
expecting light-front helicity conservation because the non-perturbative wave function is a 
mixture of different light-front-helicity states. 

If the value of becomes asymptotically large, the effects of the non-perturbative wave 
function may disappear and perturbative effects, which do respect helicity conservation, 
could take over. But for the present, we must take helicity non-conservation as a given and 
it will be worthwhile to consider the implications of this non-conservation. 



V. PROTON FORM FACTORS 

We obtain the form factors by using the wave function of Eq. ([I8| ) in Eq. (^). Our result 
can be expressed as 

F.m = J ^^"i^^^ ^('f^^Jj $nM^)$(Mo) (x^^^(p;,P^)|x^o^^(Pi,P2)) (TPsIKps)) (22) 
QkF,{Q') = 2M^ J ^HM^MMo) (Xo^'(p;,p'2))|Xo^'(Pi,P2)) (Tp^^IKps)) (23) 
The value of Mq is obtained by using Eq. (^), and 



HMo)^J^^^HMo) (24) 



El 



\ ' ^3 - Mo - E,,. (25) 



Numerical evaluations of these equations for large values of were presented in Ref . . 
Our aim here is to understand the essential features of the different Q^-dependence of Fi 
and F2. The expressions for the form factors differ only by the presence of the last factor 
(TPsIIPs)) (Tp'sITPs)! with Fi depending on the non-spin-fiip term, and F2 on the spin- 



flip term To proceed we evaluate these overlaps, using the Melosh transformations 

with = — K^, and = — + r/qj.. One computes the product of the matrices 
^Af(P3)'^if (Ps)- The upper diagonal element (non-spin-flip term) appears in the expression 
for Fi, and the upper off-diagonal element (spin- flip term proportional to cr) determines F2. 
The evaluation is simplified by realizing that integration over d'^K± causes terms linear in 
the component of K which are perpendicular to q to vanish. 

To be definite, we take q_L to lie along the x-direction, so that we find 
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(tPsltPa) 
(tPsliPa) 



[(m + (1 - r])Mo){m + (1 - r/)M^) + Kj ~ r]QK,] 
((m + (1 - r])Moy + Kl) ((m + (1 - r/)M^)2 + (K^ - r^qj 
[r7g(m + (1 - r7)Mo) + (1 - r7)(M^ - Mo)Kj 



(26) 



'((m + (1 - r/)Mo)2 + Ki) ((m + (1 - r/)M^)2 + (K^ - ryq^)2) 

We may understand the qualitative nature of the ratio using the notion that the 

value oi Q = -\/qi can be much larger than the typical momenta, of order /? = 560 MeV, 
which appear in the wave function. Then for Q ':$> j3, we may approximate Eq. (H) by 



Q^ 



1 — 1] 



and take the terms of the bracketed expressions of Eq. 
dominant. Then, using Eqs. (^) and ( [27| ) in Eqs. ( P^ . 
contains an explicit factor of Q, and 



QkF^ 



2M, 



N- 



{ri{m + (1 - r7)Mo) + ^r]{l~r])K^) 



{-K^ + (m + (1 - r/)Mo) - r^)) 



(27) 

6]) which are proportional to Q as 
we see that each of Fi and QF2 



(28) 



where the expectation value symbols abbreviate the operation of multiplying by the re- 
maining factors of Eqs. (|22| , p3D (without approximation) and performing the necessary six- 
dimensional integral. The terms r/, Mq, can be anticipated to have an expectation value 
independent of Q^, so the ratio is anticipated to be constant. 

The exact model calculation and the approximation ( P5[ ) are compared in Fig. 1. Eq. ((P5[) 
qualitatively reproduces the constant nature of the ratio and its value. Thus the constant 
nature of the ratio is understood from the properties of the Melosh transformation, which 

here embodies the relativistic effects. 
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FIG. 1. Model calculation of Eqs. (2^,23,26) solid, vs. the approximation Eq. (|2^), dashed. 
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Equation (pH]) represents a simple quick argument which gives a constant ratio. But 
this is only a rough approximation because each of Fi_2 is over predicted by about 40%. 
Numerical work shows that neglecting the terms proportional to in both the numerator 
and denominator of Eq. (pSD leads to a different approximation: 



{ri{m+{l-ri)Mo)) 



((^ + (l_^)Mo)v'r/(l-r/)) 



(29) 



which, as shown in Fig. 2 leads to an even better reproduction of the model results for Fi 
and F2. 
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FIG. 2. Model calculation of Eqs. (22,23,26) solid, vs. the approximation Eq. (pOf), dashed. 



Equation ( poD is a better approximation because the terms involving cancel against 
terms involving the difference between Mq and its approximation (^7\). Thus it seems that 
values of going up to 20 GeV^ are not large enough to allow one to completely neglect 
other terms, and therefore also not large enough to extract the asymptotic behavior. 

This feature of not reaching asymptotic values of may be understood by examining 
the dependence of the integrands of Eqs. (|22| , |23D on the value of rj. We may write 



(30) 



and determine the important regions by examining Ii^2{fliQ^)- As shown in Figs. (0,§) the 
important contributions occur for a very narrow band of values close to 1 — 77 = ^3 = 0.145. 
The sharp peaking is maintained for all of the values of considered here, and is a central 
reason for the qualitative success of the approximations ( P5|PUD . The small factor 1 — 77 
multiplies the large factor Q appearing in Eq. (0), and suppresses the dominance of the 
terms proportional to Q. The integrands peak at X3 = 0.85, a large value (compared to 0.33, 
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expected if each quark were to carry the same momentum) which indicates the presence of 
the Feynman mechanism, and a corresponding difficulty of using simple arguments to extract 
asymptotic properties of form factors. 




0.142 0.144 0.146 0.148 0.15 0.152 

FIG. 3. Important region of integration for Fi, Eq. (|30|). The curves show the derivative of Ii 
for values of = 2, 4, 6, 8, 10 GeV^, with the larger values occurring for the smaller values of Q^. 




0.142 0.144 0.146 0.148 0.15 0.152 

FIG. 4. Important region of integration for F2, Eq. (|30| ) The curves show the derivative of Ii for 
values of = 2, 4, 6, 8, 10 GeV^, with the larger values occurring for the smaller values of Q^. 
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We also find that the computed value of the ratio QHF2/F1 is remarkably independent 
of the parameters of the model. For example, Fig. (^) shows that a 10% increase in the 
value of 7, Eq. (|1^, causes about a 50% decrease in the computed values of F2, but Fig. 
shows only a 5% change in the ratio. 



0.5 




FIG. 5. The effect of varying the parameter 7 which governs the power of the falloff of the wave 
function of Eq. (0). The curves for F2 are labeled by the value of 7 = 3.5 correct model value 
(solid) or 7 = 3.9 (dashed). 

The solid curve of Fig. (^ represents our final results and predictions for values of 
that are not yet measured. The dashed curve of Fig. (^ is closer to the data for the ratio, 
but it corresponds values of F2 at high which are much smaller than the data of |^ . 
Thus the original model of Schlumpf seems to account for both Fi and F2. There is a 
small disagreement (~ 15%) with the data for the ratio, which can not be fixed simply by 
varying the parameters. This small degree of disagreement between the model and the data 
is remarkable because so many plausible effects, such as configuration mixing involving both 
quark and gluon degrees of freedom and a non-perturbative variation of the constituent 

are surely present, but these do not 



quark masses are ignored. Pion cloud effects | 
seem to be significant for values of greater than about 2 GeV^. 

We also study the behavior of the ratio QHF2/ Fi for very large values of see Fig. 7. 
The constant nature of the ratio seen in previous figures is actually the result of a broad 
maximum occurring near 10 GeV^. The ratio falls for asymptotic values of Q^, the 

displayed curve can be represented by but not as predicted using perturbative QCD, which 
would be g^^Fa/Fi ~ l/log(g2) for 80 < < eOOGeVl 
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C\2 
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GeV^ 

FIG. 6. The effect of varying the parameter 7 which governs the power of the falloff of the wave 
function of Eq. The curves for Q F2/F1 are labeled by the value of 7 = 3.5 (solid) or 7 = 3.9 

(dashed).The data for 2 < < 3.5 GeV^ are from Ref. 1., and that for 3.5 < < S.SGeV^ are 
from Ref. 2. 
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FIG. 7. High Q2 behavior of QKF2/F1 
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VI. SUMMARY AND DISCUSSION 



We have seen how a simple relativistic constituent quark model accounts for both F2 and 
QF2/F1 for values of between 2 and 5.5 GeV^. The most relevant ingredient in the model 
is its attempt to use a wave function which is Poincare invariant. In such wave functions 
the helicity conservation rule is not satisfied because the non-perturbative wave function 
is a mixture of different helicity states. This feature leads to an analytic understanding, 
embodied in Eqs. (^) and (pQ]), that QF2{Q^) and Fi(Q^) have the same variation with Q^. 
The predicted value of the ratio Q F2{Q^) / Fi{Q^) ^ 0.8 for the values of up to 20 GeV^ 
and drops slowly for larger values. 

The present data set extends to = 5.5 GeV^, and its measurement of the non- 
conservation of helicity has implications for other exclusive processes involving protons. 
Helicity conservation should be not relevant if we consider proton-proton scattering at high 
momentum transfers, up to —t = 5.5 GeV^. This means large values of various analyzing 
powers can be expected. Perhaps the most interesting mystery in proton-proton scattering 
is the large value of Ajvat observed in 90° proton scattering at s ~ 20 GeV^ [lO'Hl • This 
corresponds to — t ~ 7 — 10 GeV^. If the present measurements are extended to values of 
such as these, and if the constant nature of the ratio QF2{Q'^) / Fi{Q'^) is maintained, one 
could be able to seek an explanation of the large value of the analyzing power in terms of 
the non-perturbative proton wave function. 
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